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ABSTRACT 

We show that spiral waves in galaxy discs churn the stars and gas in a manner that 
largely preserves the overall angular momentum distribution and leads to little increase 
in random motion. Changes in the angular momenta of individual stars are typically 
as large as ~ 50% over the lifetime of the disc. The changes are concentrated around 
the corotation radius for an individual spiral wave, but since transient waves with a 
wide range of pattern speeds develop in rapid succession, the entire disc is affected. 
This behaviour has profound consequences for the metallicity gradients with radius 
in both stars and gas, since the ISM is also stirred by the same mechanism. We find 
observational support for stirring, propose a simple model for the distribution of stars 
over metallicity and age, and discuss other possible consequences. 

Key words: Galaxy: abundances - Galaxy: kinematics and dynamics - galaxies: 
evolution - galaxies: ISM - galaxies: structure - ISM: abundances 



1 INTRODUCTION 



The metallicity of the ISM increases through the life of a 
galaxy as dying stars throw newly formed heavy elements 
into the interstellar medium (ISM). It is well-established 
that the metallicity of the ISM decreases outwards in the 
Milky Way (Wilson & Rood 1994) and other galaxies (Henry 
& Worthey 1999). Since stars have the metallicity of the 
material from which they were formed, we should expect: 
(i) older stars to be more metal-poor than younger ones, 
and (ii) for a given age, the metallicity of a star to be a 
decreasing function of the radius at which it was born. 

If there were no stellar migration, the stars at a given 
radius would have metallicities that varied from zero up to 
the present-day value in the ISM at that radius and there 
would be a perfect correlation between the ages and metal- 
licities of stars. Edvardsson et al. (1993) do not find a tight 
correlation between ages and metallicities of stars in the so- 
lar neighbourhood, but find instead that stars of a given 
age have a broad spread in metallicity; e.g. stars of ages 
~ 5 Gyr have -0.7 < [Fe/H] <0.2. Only the lower bound of 
the metallicity distribution decreases with age - there are 
no young, metal-poor stars in the solar neighbourhood. 

Because of a general increase of non-circular motion 
with age (Wielen 1977; Lacey 1984; Edvardsson et al.; 
Dehnen & Binney 1998), older stars oscillate more in ra- 
dius than do younger ones, causing stars from an increas- 
ing spread of 'home radii' to appear in a sample at any 
radius. Typical radial excursions for a population of stars 
with rms radial velocity a u are A ~ \/2o- u /k, where is k is 
the epicycle frequency. For old stars in the thin disc near the 



Sun A ~ 1.3 kpc, and migration over distances of this order 
proves insufficient to explain the weakness of the correlation 
between age and metallicity found by Edvardsson et al. 

The spread in the age-metallicity relation therefore re- 
quires either that the metallicity of the ISM was much 
less homogeneous in the past than it is today, a distinctly 
unattractive possibility, or that stars have migrated by more 
than their epicycle size from their radii of birth, as proposed 
by Wielen, Fuchs & Dettbarn (1996). Since the home radius 
is determined by the star's specific angular momentum, L, 
about the Galaxy's symmetry axis, the latter possibility re- 
quires L to be changed. 

Changes in angular momentum can arise only from non- 
axisymmetric forces, and two sources have been identified: 
molecular clouds (Spitzer & Schwarzschild 1953) and spi- 
ral arms (Barbanis & Woltjer 1967), although the spiral 
supporting response of a disc to a massive object orbiting 
within it makes this distinction somewhat artificial (Julian 
& Toomre 1966). Subsequent development of these ideas 
(Lynden-Bell & Kalnajs 1972; Carlberg & Sellwood 1985; 
Binney & Lacey 1988; Toomre & Kalnajs 1991; Fuchs 2001) 
has shown that changes in the distribution of angular mo- 
mentum among the stars should always be accompanied by 
an increase in random motion, or 'heating', of the stellar dis- 
tribution. While the random velocities of stars have tended 
to increase over the lifetime of the Galactic disc, as noted 
above, even old disc stars today have A/R <C 1 and this fact 
indicates rather modest changes to the initial distribution of 
angular momentum. 

A small change in the distribution of angular momen- 
tum does not, however, forbid large changes to the angular 
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momenta of individual stars. In this paper, we show that 
spiral waves in discs do alter stellar angular momenta while 
largely preserving the overall distribution of L and causing 
little increase in random motion. We argue that churning of 
the disc by spiral waves in this manner also moves clouds 
of interstellar gas, and we discuss observational support for, 
and some consequences of, this prediction. The associated 
radial migration of stars and gas must reduce the metallic- 
ity gradient in galaxy discs and modify radically the distri- 
bution over metallicity of the stars that one now finds at a 
given radius. 

As previously noted, the spiral response of a galactic 
disc to a co-orbiting mass clump blurs the distinction be- 
tween scattering by spirals and by mass clumps. In the 
present paper, we consider the collective response of the disc 
to such mass clumps (Toomre & Kalnajs 1991) as adding 
to the transient spirals. The direct scattering of individual 
stars by clouds is important, since it probably determines 
the shape of the local velocity ellipsoid (Carlberg 1987; Jenk- 
ins & Binney 1990), but we show in the Appendix that the 
angular-momentum changes induced by star-cloud scatter- 
ing are negligibly small. 

We begin by summarizing the theory of scattering of 
stars by spiral waves. We next present simulations which 
demonstrate that angular momentum exchanges at corota- 
tion, which have previously received little attention, are in 
fact the principal driver of radial migration in galaxies. 



2 SCATTERING BY WAVES 

In this section we obtain relations between the angular 
momentum changes induced by spirals and the associated 
changes in the random velocities of stars that will be useful 
in Section 3, where we analyze numerical simulations. 

In the rotating frame of a steady spiral perturbation 
there is an 'energy' invariant, Jacobi's integral Ej, which 
is defined as [Binney & Tremaine (1987; hereafter BT) eq. 
(3-88)] 

Ej = E- n p L, (1) 

where E and L are the specific energy and z-angular mo- 
mentum of the star in a non-rotating frame and Q p is the 
pattern speed of the non-axisymmetric perturbation. Hence, 
changes in energy and angular momentum are related by 



AE 



D. P AL. 



(2) 

We now determine the fraction of the energy incre- 
ment, AE, which appears as random motion; the remainder 
changes the energy associated with circular motion. If Jr is 
any parameter that quantifies radial kinetic energy, we can 
obviously write 

dE=^djR + §dL. (3) 
oJr oL 

If Jr is chosen to be the 'radial action', the partial deriva- 
tives in this expression become the angular frequencies, lor 
and f2, of a star's radial and azimuthal motion [BT eq. (3- 
150a)]. In the epicycle approximation, lor = n the classical 
epicycle frequency, Jr is simply Er/k, and the radial kinetic 
energy Er is related to the maximum radial speed around 



an epicycle by Er = ±(v fl , m ax) 2 . 
equations ^ and (^), we obtain 



Eliminating AE between 



Eccentric orbits Circular orbits 




Figure 1. Classical Lindblad diagram. See text for a description. 
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(4) 



While changes in L at corotation (Q — Q p ) do not cause 
changes in Jr, those away from corotation do. 

The classical Lindblad diagram, Fig. [jj illustrates the 
physical origin of this relation. Possible orbits occupy the 
upper left half of this diagram, and are bounded by the locus 
of the circular orbit of angular momentum L. Equation 
states that a steady spiral wave moves stars on lines with 
slope £l p . Since the tangent to the locus of circular orbits 
at corotation has this slope, Jr, which is a measure of the 
distance from the circular-orbit curve, does not change to 
first order when a star corotating with the wave is scattered. 
Scattering at more general locations, does result in energy 
being exchanged between random and orbital motion, and 
stars on near circular orbits must be preferentially scattered 
into non-circular orbits. 

When a star is in Lindblad resonance with an m-armed 
spiral, its frequencies satisfy [e.g. BT eq. (3-122)] 



lur = ±m(fi p — H), 



(5) 



where the upper sign applies at the outer Lindblad reso- 
nance (OLR) and the lower sign applies at the inner res- 
onance (ILR). Combining equations ([|) and (^) we obtain 
the pleasingly simple result 



AJr = =p — AL (Lindblad resonance). (6) 
m 

Here the upper sign applies if the star lies at ILR (fi > f2 p ) 
and the lower sign applies at OLR. This relation is exact 
for nearly circular orbits, and remains the dominant term 
for moderately eccentric orbits where higher resonances also 
contribute (Lynden-Bell & Kalnajs 1972). 



3 SINGLE SPIRAL 

As a test of the spiral wave theory presented in §2, we begin 
by studying what happens in a simulation having a single 
transient spiral wave. We later describe more realistic sim- 
ulations with a succession of spiral waves. 
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Figure 2. The later part of the growth and subsequent decay of an isolated spiral mode in a disc that was seeded with a groove. The 
radius of the circles is 15i?i and only one particle in 120 is plotted. 



3.1 Choice of model 

We wish to construct a simulation that supports a single spi- 
ral wave standing out clearly from the noise; to achieve this 
goal, we follow the procedure described by Sellwood & Kahn 
(1991, §3.3). Mestel's (1963) disc, in which the circular speed 
is everywhere Vb, suffers only from global one-armed insta- 
bilities when the disc is almost fully self-gravitating (Zang 
1976; Toomre 1977, 1981; Evans & Read 1998), and becomes 
stable to all global modes when only half the disc mass is 
active (Toomre 1981). These remarkable stability proper- 
ties have been confirmed in JV-body simulations (Sellwood 
& Kahn 1991; Sellwood & Evans 2001). An appropriate dis- 
tribution function (DF) for this disc is (Zang 1976; Toomre 
1977) 

CL , e l-E( q +l)] 



Table 1. Numerical parameters used in the two simulations 



fo{E,L) = -, 



(7) 



2[l + (L i /L)4][l + (L/L )«]- 

The constant C [given in BT (eq 4-162)] yields the correct 
surface density when fo is integrated over all velocities with 
the factors in the denominator set to unity, while the pre- 
factor halves the active mass. We set the parameter q = 
11.44, which in an infinite half-mass disc would yield Q = 
1.5. The radial extent of the otherwise infinite disc is limited 
by the inner and outer tapers, represented by the factors 
in the denominator; these tapers are gentle enough not to 
provoke instabilities on their own (Toomre 1981). The value 
of Q is raised where the tapers kick in, since they reduce the 
active mass density by more than the velocity dispersion. 
As the tapers also introduce length scales into the disc, we 
adopt the central radius of the inner taper, Ri — Li/Vo, as 
our length unit, and set L Q = 15Li. We also adopt Vb as our 
velocity unit, so Ri/Vo is our unit of time and Li our unit 
of angular momentum. 

Sellwood & Kahn (1991) showed that the introduction 
of a narrow axisymmetric feature into the equilibrium DF 





Model S 


Model U 


Grid size 


230 x 256 


110 x 128 


Softening length 


0.05Ri 


O.lRi 


Active components 


m = 2 only 


< m < 4 


Particle number 


6 X 10 5 


10 6 


Time step 


0.005i?i/Vb 


0.005fli/Vb 



destabilizes the disc in a predictable manner. We therefore 
multiply /o by a function of L and employ the distribution 
function 



f(E,L) = f (E,L) 



1 + 



{L-L«y + wi 



(8) 



We chose the parameters of the Lorentzian function to be: 
L, = 6.5Li, P = —0.4, idl = 0.2Li. While this represents 
a groove in the angular momentum density, it causes only 
a mild depression in the surface density over a broad range 
of radii because the particles move on epicycles; the mean 
epicycle radius A ~ 0.28.R. 



3.2 Numerical procedure 

We set up the model with considerable care to ensure an ac- 
curate initial equilibrium, and to keep random fluctuations 
to a minimum. We adopt the quiet start procedure described 
by Sellwood & Athanassoula (1986) to select the initial co- 
ordinates and place the particles at equal angular intervals 
around rings. 

We use a highly-efficient PM method (Sellwood 1981, 
1983) to compute the inter-particle gravitational forces, with 
a small softening parameter. Particles move freely over a 
two-dimensional polar lattice of {Nr,N^) points, which is 
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AE/AL 

Figure 3. The distribution of the ratio AE/AL of the changes in 
integrals up to t = 300 in the disc that was seeded with a groove. 
The vertical dashed line marks the measured pattern speed of 
the spiral mode. The shaded histogram shows only those particles 
with initial L > 2, while the unshaded shows all particles. 

used to tabulate the gravitational field. The grid points are 
spaced at R = O.liii [exp(an) — 1], where a = 2n/N < j > and 
< n < Nr. The mass of each particle is spread over the 
four nearest grid points and we employ a standard Plum- 
mer softening formula to further smooth the forces at short 
range. Radial and azimuthal acceleration components are 
computed directly and tabulated, avoiding the need to differ- 
ence a potential function, and we use bi-linear interpolation 
to evaluate the acceleration components on each particle. 
Our choices of numerical parameters in this run, model S, 
are given in Table [[]. 

Forces from the particles are evaluated every time step. 
Particles in the radial range 0.5 < R < 1 are integrated using 
the given time step, while beyond R = 1 time steps increase 
by factors of 2 for every factor of 2 increase in radius, with 
forces arising from particles in outer zones being determined 
by interpolation in time when needed to advance the motion 
of those on shorter steps - see Sellwood (1985). Time steps 
for the few particles inside R = 0.5 are sub-divided by fac- 
tors of 2 as they approach the centre, but the force field in 
this region, which is anyway dominated by the fixed central 
attraction, is not updated any more frequently. 

The groove will excite a separate instability for each 
azimuthal wavenumber, m (Sellwood & Kahn 1991). We 
therefore restrict the disturbance forces from the particles 
to those arising from the m — 2 component of the density 
distribution. The axisymmetric part of the field is held fixed 
throughout. 




0.0 



4 8 12 

R 

Figure 5. Dispersion in radial velocity as a function of radius in 
the disc that was seeded with a groove at t = (light) and t = 300 
(dark). The vertical lines mark the positions of corotation (solid) 
and the Lindblad resonances (dotted) for near-circular orbits. 

shows this quantity evaluated for only those particles with 
initial L > 2Li; a heavy tail to large values appears only 
when particles with small L are included (unshaded) , show- 
ing that the tail is associated with the formation of the inner 
bar. 

Fig. ^| shows the mean and the 20th and 80th percentiles 
of the distribution of AL values as a function of the initial 
L; corotation and the Lindblad resonances for the spiral are 
marked. The left-hand panel shows the changes for all par- 
ticles, the right-hand panel shows those for the 20% of par- 
ticles with the smallest epicycle energies. The changes are 
large, especially near corotation. AL increases with distance 
from corotation at such a rate that particles are transferred 
to the other side of corotation: the line of slope —2, on which 
particles move symmetrically across corotation, is marked. 
Thus the net effect of the large individual changes in L is to 
cause particles to change places. 

Changes away from corotation are clearly much smaller. 
Some are associated with the Lindblad resonances, but those 
near the centre are perhaps somewhat confused by the de- 
velopment of the bar. 

The theory in Section 2 asserts that the heating associ- 
ated with a given change AL tends to zero as we approach 
corotation from either side. Fig. ^confirms this prediction by 
showing the dispersion in the radial velocities as a function 
of radius at the initial moment and at t = 300. It is evident 
that the large changes in L at corotation produced no sig- 
nificant heating whereas the comparatively small changes at 
the inner Lindblad resonance caused considerable heating. 



3.3 Effect of one spiral disturbance 

Fig. ^ shows the later evolution of a single growing spiral 
mode which saturates at about time 270 and then decays 
leaving a weak bar-like feature. Fig. ^ shows the distribu- 
tion of values of the ratio AE/AL, which eq. (^) implies 
will equal Q p . The sharp peak in this distribution coin- 
cides extremely well with the measured pattern speed of 
the mode, Q p = 0.141, implying that the vast majority of 
changes are induced by this spiral. The shaded histogram 



3.4 Orbit trapping at corotation 

The physics of the changes at corotation is easy to under- 
stand. We show, in Fig. ^, orbits in a Mestel disc when 
a steady two-armed spiral perturbation is imposed. Parti- 
cles initially have velocity of magnitude Vo in the tangential 
direction and are followed until their orbits approximately 
close on themselves. Those inside corotation overtake the 
wave; they gain angular momentum as they fall into the spi- 
ral arm which causes them to move to larger radii and slow 
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Figure 4. Left panel: the full curve shows the mean change in L to t = 300 for all particles of the given initial L. The shaded area is 
bounded by the 20th and 80th percentiles of the particles when ordered by magnitude of AT. The dashed line has slope —2, to show the 
locus of particles whose changes would be symmetric about corotation. The vertical lines mark the positions of corotation (solid) and 
the Lindblad resonances (dotted) for near-circular orbits. 



their drift relative to the wave, enabling large changes in 
L to develop. As the figure shows, particles initially close 
enough to the radius of corotation, interact so strongly that 
their speed relative to the wave is reversed. These are called 
'horseshoe orbits' (Goldreich & Tremaine 1982). A parti- 
cle that has moved out across corotation in this way subse- 
quently slips backwards relative to the wave and eventually 
falls backwards into the other arm. As it falls, it loses angu- 
lar momentum, moves to smaller radii and eventually pulls 
ahead of the wave. These successive episodes of angular- 
momentum gain and loss cause particles on horseshoe or- 
bits to move around the maximum of the spiral potential as 
shown. 

Section 3.3.3(b) of BT describes an analytic approxi- 
mation to such motion in a sinusoidally varying potential of 
time-independent amplitude *o- In this approximation, the 
phase variable ip = 2{4> — m ax), where max is the azimuth 
at which the potential peaks, obeys the pendulum equation 



ip = -i 
where 
2 

P= R 



' suit/;, 



-B 



With A and B the Oort constants; 
A = -B = n/2. Eq. (|) admits 



(9) 



(10) 

for the Mestel disc, 
an energy invariant 



E p — r;ip 2 — p 2 cosip. Orbits with E p > p 2 circulate, while 
the horseshoe orbits have Ep < p 2 and librate. The periods 
of the horseshoe orbits increase with E p from the minimum 
value [BT eq. (3-127b)] 




Figure 6. Orbits in a spiral perturbation. The background model 
is Mestel's disc and the perturbing potential is <E>i = "To cos[2(0 — 
00 - &pt)] where <t>o(R) = 2(R/R - 1) + tt/2. The locus of the 
deepest part of the wave's potential is marked by dots. 



Tmin — TtRo 



-B 



A\t>n 



(11) 



that is associated with the orbits that make the smallest 
excursions in R. For the orbit on which R varies most, the 
changes in R and L are given by [BT eq. (3-129)] 



AT? = 2 



-AB 
AL = 2ARAR. 



(12) 



The spiral perturbation in the simulation is transient, 
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Figure 7. Part of the time variation of the amplitude of the 
m = 2 component of the gravitational potential, In&oI, measured 
at R = 7.1-Ri, the radius of corotation of the spiral shown in 
Fig. § 

rather than steady. Clearly, scattering from one side of 
corotation to the other cannot be completed if the non- 
axisymmetric potential is at full strength for less long than 
0.5T m i n . On the other hand, if the disturbance persists for a 
time much greater than 0.5T m j n , many particles will be scat- 
tered both across and back, leading to smaller net changes. 

The time variation of Wo near the radius of corotation 
is shown in Figure |?]. From t = to t ~ 250, the disturbance 
potential grows exponentially, with growth rate 7 = 0.047. 
At t ~ 250 growth slows and from t ~ 270 the disturbance 
potential fades even more quickly than it grew. Rough es- 
timates from Fig. [7] are that the peak amplitude of the 
perturbing potential Wo — 0.35 and its duration T on ~ 20 
time units. With eq. (|ll]) these values yield T m i n ~ 38, so 
T on ~ 0.53T m i n . Hence, as expected, particles cross the res- 
onance at most once. Inserting our value for Wo in eq. ( |l2"| ) 
we find that the largest change in angular momentum is pre- 
dicted to be AL — AR ~ 2.4, in excellent agreement with 
the observed magnitude of the changes around corotation in 
Fig.g 

This picture of the instability's physics can only be ap- 
proximate since it is based on orbits in a steady rather than a 
time- varying potential. None the less, we have shown that it 
yields a semi-quantitative understanding of the simulations. 

Moreover, it predicts that particles on nearly circular 
orbits will experience larger changes in L than particles on 
more eccentric orbits: the angular velocity with which the 
latter circulate about the galactic centre varies significantly 
as the particles oscillate radially, so it is impossible for such 
particles to hold station with respect to a steadily rotating 
wave. The right-hand panel of Fig. ^ confirms this prediction 
by showing the mean and spread of changes in L for the 20% 
of particles with the smallest epicycle energies. One sees that 
the changes are ~ 50% larger than those for the disc as a 
whole. 

The emergence of horseshoe behaviour may be the prin- 
cipal change that limits the amplitude of a spiral mode. 
Standard linear stability theory (e.g. Kalnajs 1971, 1977) as- 
sumes an infinitesimal disturbance and computes the orbital 
response as small departures from the unperturbed orbits. 
The linear mode is the self-consistent disturbance for which 



the sum of the orbit responses is equal to that required to 
generate the disturbance potential with a particular pattern 
speed, fij,, and growth rate, 7. At any finite amplitude there 
will be particles on horseshoe orbits, but the existence of 
these orbits will not have a significant impact at first, be- 
cause (a) they will be small in number, and (b) they will have 
long periods, and first-order perturbation theory will give a 
good account of the motion of a particle on a horseshoe orbit 
until the particle reaches the point at which tp = 0. As the 
mode grows the period of any trapped orbits decreases, and 
the width of the horseshoe region grows exponentially; thus 
the number of particles that develop horseshoe behaviour 
grows exponentially and the time for which linear perturba- 
tion theory remains valid for them becomes short. When the 
condition 7T m i n > 1 is first violated, particles on horseshoe 
orbits that have passed through tp = are suddenly impor- 
tant contributors to the overall dynamics, and we have to 
expect linear perturbation theory to fail. Thus many parti- 
cles reverse their direction of motion in the rotating frame 
rather abruptly, which leads to a rapid decrease in the dis- 
turbance potential as the density enhancement disperses on 
a timescale >0.5T m i n (~ 20). These ideas are consistent 
with the behaviour shown in Fig. 

3.5 Comparison with other work 

Lynden-Bell & Kalnajs (1972; hereafter LBK) presented 
an analysis of the angular-momentum exchanges between 
a growing wave and a disc that went far beyond the ele- 
mentary considerations of Section 2. Their widely quoted 
result is that in the limit of vanishing growth rate angular- 
momentum changes occur only at resonances. Carlberg & 
Sellwood (1985; hereafter CS), focusing on heating, extended 
their analysis to transient waves. 

LBK averaged the change in L over all initial phases 
of the star. Clearly, orbit-averaging makes it impossible to 
discover how wide is the spread in AL at a given initial 
value of L. However, one might hope that it would correctly 
predict the net change in L, and qualitatively it does; their 
results for growing waves correctly predict that resonances 
are broadened, and that stars just inside corotation should, 
on average, gain L while there should be a net loss just 
outside. 

Since the second-order perturbation theory developed 
by LBK and CS assumes that changes to all quantities are 
small, their analyses implicitly exclude the trapping pro- 
cess associated with horseshoe orbits (Fig. fy, but may give 
an adequate description of the behaviour outside the trap- 
ping region. The signs of the angular momentum changes in 
Fig. ^ are in agreement with LBK and CS, but the largest 
changes in L are associated with order unity changes to an 
angle variable and are therefore outside the regime where 
their predictions hold. Furthermore, any perturbation anal- 
ysis which seeks to expand changes in the dynamical vari- 
ables in powers of the perturbing potential Wo, will have 
difficulty in recovering the result AL ~ v/jWoj. Horseshoe 
orbits become important in the simulation both because the 
perturbation is strong and because it is of short duration. 
A stong perturbation causes the horseshoe region to have a 
significant width, while a more enduring perturbation would 
allow particles to recross corotation and end up with essen- 
tially unchanged angular momenta. 
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Figure 8. Part of the evolution of the unconstrained simulation. The spiral activity is very mild at all times, the pattern visible at time 
1170 is perhaps the strongest of all. The radius of the circles is 21R; and only 1 particle in 200 is plotted. 



In an influential paper Wielen (1977) showed that the 
increase in the velocity dispersion of K and M stars with age 
could be reproduced by a model in which stars diffused in 
velocity space with diffusion coefficients that are isotropic 
and constant. In reality, even in the absence of scattering, 
the velocity of a star changes radically in the course of a half 
an epicycle period, so a picture based on velocity space is of 
limited value. When stars are viewed as diffusing in integral 
space (Binney & Lacey 1988), one finds both that diffusion 
is extremely anisotropic, and that the diffusion coefficients 
vary from point to point in the space in a manner which 
contradicts Wielen's simple picture. Indeed, star-wave scat- 
tering drives diffusion that is as anisotropic as it logically 
can be, in that stars diffuse along lines in two-dimensional 
integral space (Section 2). Neglecting such effects leads to 
seriously mistaken results. 

It seems possible that the horseshoe behaviour we have 
identified is the effect that Zhang (1996) describes as a "col- 
lisionless shock" as stars pass through spiral arms. She also 
claims additional scattering in this event, but any scattering 
at corotation must occur without heating the disc. 



4 UNCONSTRAINED SIMULATION 

We next describe a more realistic simulation which displays 
a succession of very mild spiral waves. We again use a half- 
mass Mestel disc with initial Q = 1.5 (in the absence of 
tapers). This model has no initial groove and we place our 
one million particles at random azimuths, so that the model 
starts with shot noise. We also include forces from azimuthal 
wavenumbers < m < 4, so that there are often significant 
patterns with different rotational symmetries present at the 
same time. The numerical details, summarized under Model 
U in Table [j], differ in insignificant respects from those in 
Model S described above. 

The simulation, shown in Fig. pi displays a succession 




4 8 12 16 20 

Radius 



Figure 9. The light shaded region shows the initial radial velocity 
dispersion in model U, the dark shaded region shows the change 
by t = 1200. The two curves and right-hand scale show the local 
stability parameter Q at the same times. 



of very mild spiral waves that gradually heat the inner disc, 
while the outer disc has not heated much by t — 1200. The 
change in the radial velocity dispersion between the start 
and the end of the simulation is illustrated in Fig. ^; it more 
than doubles for stars in the inner disc. The stability pa- 
rameter, Q, rises to Q > 2 everywhere by t = 1200, when 
spiral activity has diminished but not completely ceased. 

Fig. [it] shows as a function of radius the power in the 
Fourier transform with respect to time of the m = 2 com- 
ponent of the density distribution. The power is strongly 
concentrated around a number of discrete frequencies u>i. 
The Lindblad resonances for each frequency occur where 
the horizontal ridge intersects the dashed curves and the in- 
tersection with the full curve marks corotation. It is evident 
that there is always significant power extending from the 
ILR to CR for each wave, and the power sometimes extends 
out to the OLR. 
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Figure 10. The power spectrum of the m = 2 component of the 
density variations in the unconstrained simulation. The dashed 
curves mark the loci of ILR and OLR for the given frequency, 
while the full curve shows the locus of CR. 




L(0) 

Figure 11. The change in L to t = 1200 plotted against the 
initial L for one particle in 11 of those in the simulation with 
recurrent transient spiral structure. 

Fig. [ll] shows for each particle the difference, AL, be- 
tween the initial and final values of L for that particle plotted 
against the initial value of L. Changes in L of 50 percent are 
not uncommon. Fig. |t| which is an unsharp-masked version 
of Fig. [ll], reveals small-scale variations in particle density 
that in Fig. [ll] are obscured by the overall decline in the par- 
ticle density from the line AL = 0. Several diagonal ridges 
of enhanced particle density are apparent. Each ridge is pre- 
sumably the signature of a spiral feature of some frequency 
tUi, and is equivalent to the shaded strip in Fig. ^ as it runs 
from upper left to lower right through CR. Thus each di- 
agonal ridge in Fig. [l2| is created by particles on horseshoe 
orbits crossing corotation for some Ui. 

Figure [13 shows the extent of mixing within the disc 
that spiral activity causes by showing several distributions 



3 r 




L 

Figure 12. An unsharp-masked version of Fig. [ll] reveals ridges 
in the particle distribution that extend diagonally from top left 
to bottom right. 

in final home radius of particles that initially have essentially 
identical home radii. The distributions are wide: particles at 
intermediate radii can be moved by a factor 2 or more in ei- 
ther direction. The initial mean epicycle size, A ~ 0.28iihomc 
throughout most of the disc, doubles for stars in the in- 
ner disc, while those at large radii are scarcely changed by 
t = 1200 (Fig. §). 



5 DISCUSSION 

5.1 Stellar migration 

Our simulations suggest that stars can migrate significant 
distances within discs, whilst remaining on nearly circular 
orbits. To decide how important this effect is for real galax- 
ies, we need to scale the numerical results to take into ac- 
count that galaxies are typically three times older than the 
duration of our simulation, and to allow for possible differ- 
ences in the amplitude |^o| of the non-axisymmetric com- 
ponents of the potentials. In Section 3.4, we argued that the 
size of a typical step in radius is proportional to ^/|^o|, so, 
if successive steps are uncorrelated and on average equally 
sp aced in time, the net distance migrated should scale like 
y/ t\$fo\, where t is the age of the disc. 

Near-infrared photometry provides the obvious way to 
estimate \I>o- Rix & Rieke (1993) have argued that the con- 
tribution by massive stars to K-band luminosity is small 
and patchy, so the A"-band light is a reliable measure of 
the stellar mass distribution in galaxies on large scales. Rix 

6 Zaritsky (1995), Gnedin et al. (1995), Block & Puerari 
(1999) find that in the A'-band the arm-interarm surface- 
brightness contrast is generally >30%. For comparison in 
our unconstrained simulation the arm-interarm contrast in 
mass density is at most ~ 15%. The extent of radial mi- 
gration must also be a sensitive function of the disc-mass 
fraction, because both the arm-interarm contrast and the 
magnitude of the potential fluctuation that is produced by 
a given contrast increase with disc mass- fraction. A variety 
of arguments indicate that in high surface-brightness disc 
galaxies, more than half of the mass interior to two scale 
lengths lies in the disc (Debattista & Sellwood 1998, 2000; 
Weiner, Sellwood & Williams 2000; Kranz, Slyz & Rix 2002; 
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Figure 14. Oxygen abundances of 31 HII regions from Shaver et 
al. (1983; squares) and from Vilchez & Esteban (1996; triangles). 
The original data have been rescaled to Ho = 8kpc. The Orion 
nebula, from Shaver et al., is marked by a cross, and the Sun, with 
O abundance 8.93 from Anders & Grevesse (1989), is marked by 
a circle. The line is the linear least-squares fit to the data for 
R > 5.5 kpc. 
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Figure 13. The distributions of home radii at t = 1200 of par- 
ticles starting with initial home radii within the range indicated 
by the vertical dotted lines. Typical epicycle radii are initially 
~ 0.28i?home; t ne y do not increase much in the outer disc but 
rise in the inner disc to ~ ().5/<i u)Im , by t — 1200. 



?r 



-0.5 



Binney & Evans 2001). In our simulations exactly half of the 
mass resides in the disc, so we are probably underestimating 
the extent of radial migration by a factor of a few. Widen- 
ing the distribution shown in the left-hand middle panel of 
Fig. |l^ by even a factor 1.5, one is led to the conclusion that 
old stars formed in the solar neighbourhood should be scat- 
tered nearly uniformly within the annulus from R = 4 kpc 
to R= 12 kpc. 

If the radial migrations of individual stars were equiv- 
alent to the random steps taken by particles in a diffusive 
medium, the galactic disc would spread quite rapidly be- 
cause the radial steps are rather large. In the present case 
however, any spreading of the disc will actually be rather 
slow because the large migrations that occur each side of 
corotation almost exactly cancel by conservation of angu- 
lar momentum; only the relatively small transfer of angu- 
lar momentum outwards from the ILR to corotation en- 
ables the disc to spread radially. The inability of angular- 
momentum exchanges at corotation to induce radial spread- 
ing is reflected in their inability to heat the disc significantly 
[eq. (Q)], since when a disc spreads, random energy is in- 
evitably released (LBK; Lynden-Bell & Pringle, 1974). 



5.2 Abundance distributions 

The really important effects of radial migration lie in the 
field of chemical evolution. Fig. [l4| plots measured values 
of O/H against Galactocentric radius for HII regions from 



10 



T/Gyr 



Figure 15. The distribution of a sample of 189 nearby stars in 
metallicity and age. The data are from Edvardsson et al. (1993). 
The circle marks the position of the Sun. 



Shaver et al. (1983) and Vilchez & Estaban (1996) rescaled 
to Ro ~ 8 kpc. The line is the linear least-squares fit to the 
data for R > 5.5 kpc. The dispersion in [O/H] about the 
mean line is 0.16 dex, which is not significantly larger than 
the quoted errors in the data of Shaver et al. and Vilchez 
& Estaban. Hence, it is widely believed that the intrinsic 
dispersion in [O/H] at fixed 7? is < 0.1 dex (Edmunds 1998) 
and at each radius in the disc we may suppose that the ISM 
has a well-defined metallicity. 

The squares in Fig. |l^ show the distribution in [Fe/H] 
and age of 189 stars from the sample of Edvardsson et al. 
(1993). The Sun is plotted as a circle. Since Edvardsson 
et al. selected stars to obtain wide ranges in r and [Fe/H], 
the distribution of points has significant selection bias in the 
sense that the spread in age and [Fe/H] is unrepresentatively 
broad (see §3.2 of Edmunds 1998). None the less, the figure 
shows that age and metallicity are weakly correlated if at 
all; the only clear effect is that there is a dearth of young, 
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r/Gyr 

Figure 16. The density of points in the (log(Z), log r) plane pre- 
dicted by the models of chemical evolution and stellar migration 
described in the text. The curve shows the local metallicity of the 
ISM according to equations ( |l3| ) and 



metal-poor stars, presumably because the ISM at all star- 
forming radii has been for some time more metal-rich than 

[Fe/H] 0.4. 

Radial migration enables old stars that formed from 
promptly-enriched gas at small galactocentric radii to ap- 
pear in a solar-neighbourhood sample, thus weakening any 
correlation between age and metallicity. Migration from out- 
side Ro also allows a smaller number of young metal-poor 
stars to be in such a sample. A simple quantitative model 
illustrates these phenomena. In this simple model, at time 
t and radius 7? stars form at a constant rate from an ISM 
having a metallicity rising linearly with time as 



Z(R,t) = Z G + -[Z(R,t ) - Z G ], 
to 



(13) 



where Zq is the metallicity of protodisc material and to = 
12 Gyr is the age of the disc. At the present time the metal- 
licity of the gas is given by the line in Fig. 14 



log[Z(R,t )} 



(I- 1 ) 



where 



a = -0.17 
6 = 0.96. 



(14) 



The density of stars in the (log Z, log r) plane is then given 
by 

dR 



N Ro (t, Z) oc TZRe~ R/Rd P T (Ro\R) 



dZ 



(15) 



Here the surface density of newly formed stars is assumed 
to be proportional to e~ H ^ Hd with R d = 0.28i?o (Drimmel 
& Spergel, 2001), while P T (R \R) is the probability that 
a star formed at radius R migrates to -Ro in time r, and 
R(to — t, Z) is the radius at which the ISM had metallicity 
Z at time to — r. We calculate illustrative predictions for the 
case in which P T is 



P T (Ro,R) = (27va 2 T 
with 
o> = Ro 



-1/2 



exp 



(Ro - R) 2 

2cr? 



(0.16) 2 + (0.4) 2 - 
to. 



1/2 



(16) 



(17) 



The first term in the square brackets represents the radial 
migration associated with epicyclic motion for a typical thin- 
disc radial velocity a u = 34 km s _1 (Table 10.4 of Binney & 



Merrifield 1998), which we assume for simplicity to be time- 
independent. Fig. [Hi] shows the predicted distribution that 
one obtains from these formulae with Zg = O.IZq, which 
is the approximate lower boundary on the populated zone 
in Fig. ha. The distribution of Fig. [16| h as features in com- 
mon with the observational plot, Fig7Tl5l In particular both 
plots show a dearth of stars around (2, —0.5), a high popula- 
tion density around (8, 0) and an approximately horizontal 
upper boundary on the populated zone. The observational 
plot does not show the same extremely high density of stars 
around (10, —0.75) that is evident in the predicted distribu- 
tion. This density peak is large due to our assumption that 
at early times the ISM everywhere had metallicity O.IZq, 
but it is exacerbated by the use of a logarithmic age coordi- 
nate. In the observational plot errors in age estimates, which 
are obviously largest for the oldest stars, scatter stars from 
the furthest realistic ages to unrealistically high ages, thus 
reducing the density of stars around r = to- 

Wielen et al. (1996, hereafter WFD) argue that the Sun 
is a beautiful example of stellar migration. It is on an unusu- 
ally circular orbit for a star of its age (4.5 Gyr), with home 
radius about 200 pc outside Ro and pericentre about 140 pc 
inside Ro. WFD note that the Sun is not only more metal- 
rich than the average nearby star of comparable age by 
0.17 dex (Fig. ||), but is actually more metal-rich than the 
present-day local ISM (as represented by the Orion molecu- 
lar cloud, which is marked by a cross in Fig. |l4| ). They argue 
from these facts that the Sun was born R ~ (Ro — 2kpc). 

Fig. [l6| clearly shows that migration brings into the so- 
lar neighbourhood many more stars that are more metal rich 
than the local ISM than it imports more metal-poor stars: 
metal-rich immigrants are imported by arms with corotation 
inside Ro that effect exchanges between the solar neighbour- 
hood and an annulus around R ~ 5kpc, while metal-poor 
immigrants are imported by arms that have corotation out- 
side Ro and move stars from R ~ llkpc. Since there are 
many fewer stars at R ~ 11 kpc than at R ~ 5kpc, metal- 
rich immigrants outnumber metal-poor ones by a significant 
factor. The clear implication is that at R ~ 5 kpc there 
should be many stars with sub-solar metallicities. As most 
migration is accomplished with little change to the random 
motion of the stars, it is not surprising that the Sun has a 
nearly circular orbit. 

5.3 Interstellar gas 

Gas clouds are deflected by the spiral-arm potential in much 
the same way as are stars. The deflections generate in-going 
and out-going streams that occur at different azimuths and 
do not cross until large changes in radius have been achieved 
(Sellwood & Preto 2002). The resulting flow in the coro- 
tation annulus has a characteristic anticyclonic form. Our 
model predicts that the entire flow lasts no more than half its 
turnover time and will involve velocities \v r \ < | M/o | / C — BRo). 
Fridman, et al. (2001) have already called attention to flows 
of this form in NGC 3631, and derive peak radial velocities 
approaching half the orbital speed at their estimated coro- 
tation radius, again suggesting a very strong spiral pattern. 

Fig. |i| showed that stars with lower random velocities 
experienced larger changes in L than did hotter stars. Since 
gas clouds have smaller random velocities than any popula- 
tion of stars, they should experience still larger changes in L. 
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Hence, it is possible that the Orion nebula lies so far below 
the mean line in Fig. ^ because it has immigrated to the 
solar neighbourhood from R ~ lOkpc, where its metallicity 
would be typical. If the metallicity distribution of clouds at 
a given radius has non-negligible width, the theoretical stel- 
lar metallicity distribution of Fig. which is based on the 
assumption of negligible width, will be too narrow at small 
ages. Comparison of Figs ^| and [Hs] suggest that this may 
be true, but one has to bear in mind that the width of the 
observational distribution is exaggerated by both selection 
bias and errors in the measured stellar metallicities. 

5.4 Dust grains 

The impact of migration on the metallicity of interstellar 
gas must be to some extent masked by mixing following col- 
lisions between clouds of different metallicity, to produce a 
new cloud of intermediate metallicity. The chemical com- 
position of an individual interstellar grain cannot change, 
however, so the compositions of the interstellar grains that 
are currently located at a given radius should show similar 
evidence for migration to that evident in Fig. |l5|. Clayton 
(1997) presents evidence for just such an effect. 

As the disc ages, the isotopes 29 Si and 30 Si, which are 
secondary, should become more abundant relative to 28 Si, 
which is primary. Dust grains that were formed before the 
solar system have been identified in meteorites, and mea- 
surements of the relative abundances of the three Si isotopes 
(Hoppe et al. 1994a,b) confirm the theoretical prediction 
(Clayton 1988; Timmes & Clayton 1996) that the abun- 
dance ratios n( 29 Si)/n( 28 Si) and n( 30 Si)/n( 28 Si) should be 
tightly correlated in a positive sense. Three problems mar 
the beauty of the observed correlation: the slope of the corre- 
lation differs from that expected; the Sun's ratios are barely 
compatible with the correlation; and the measured ratios 
indicate that the majority of presolar grains formed from 
material that was substantially more heavily processed, nu- 
cleosynthetically, than that for which the Sun formed. Clay- 
ton (1997) points out that this last finding is similar to the 
metal-richness of the Sun relative to the local ISM because 
it implies that material that left the ISM at some time prior 
to the formation of the Sun, somehow experienced more nu- 
cleosynthesis than material that stayed in the ISM, and con- 
tinued to be enriched, right up to the formation of the Sun. 
Clayton further argues that the large changes in home ra- 
dius proposed by WFD could explain the Si isotope ratios in 
presolar grains by making it credible that they formed from 
material that had been enriched by stars that themselves 
formed ~ 2kpc interior to the birthplace of the Sun. In the 
light of our results we suggest that these grains themselves 
migrated from R ~ 4kpc to R ~ 6kpc, where they were 
incorporated into the protosolar cloud. 

5.5 Galactic dynamos 

Radial mixing in the ISM may also help with the well-known 
problem posed by the large-scale (ordered) component of B- 
fields in galaxies (e.g. Rees 1994) . Standard afi-dynamo the- 
ory (Parker 1955) is thought to yield too low a growth rate to 
achieve the present-day observed field strengths (Beck, et al. 
1996) from the likely seed fields. The growth-rate is propor- 
tional to the geometric mean of the rates of galactic shear 



(the Q term) and cyclonic circulation (the a term) (Kul- 
srud 1999). Current estimates of the a-term are based on 
supernovae-driven turbulence (Ferriere 1998), but the large- 
scale radial mixing discussed here may enhance the a-effect 
substantially, and thereby increase the growth-rate obtain- 
able from the dynamo. 



6 CONCLUSIONS 

We have shown that the dominant effect of spiral waves is to 
churn the disc. Stars just inside corotation swap places with 
those outside because each group surfs on opposite sides of 
the wave. We present a simplified model of this behaviour 
based on the theory of horseshoe orbits in a steady poten- 
tial. We find that transient spiral arms are at peak strength 
for long enough to produce generally only a single crossing 
of corotation by a large fraction of the stars over a broad 
swath of the disc around this resonance. The length of the 
largest steps taken in R scales as •y/j^oj, where is the 
peak amplitude of the perturbing potential. Scaling to real 
galaxies requires calibration of the near IR photometry and 
estimates of the disc mass fraction, but step sizes of 2 - 3 
kpc are expected. 

Groups of stars move in both directions, and gener- 
ally just exchange places. This exchange occurs without any 
significant radial spreading of the disc or increase in non- 
circular motions. In this respect exchanges at corotation 
differ significantly from those at the Lindblad resonances, 
which do engender heating and spreading. We emphasize, 
therefore, that the dominant scattering process in the disc 
does not alter the surface density profile of the disc. 

We argue that the peak amplitude and duration of a spi- 
ral wave are limited by horseshoe behaviour at corotation. 
The density enhancement disperses rapidly when a substan- 
tial fraction of particles on horseshoe orbits have reversed 
their motion in the rotating frame. 

The induced radial migration of stars will largely erase 
the correlation between metallicity and stellar age that is 
a clear prediction of standard chemical-evolution theory. It 
will also introduce into the solar neighbourhood large num- 
bers of stars that are substantially more metal-rich than the 
local ISM was at the time of their birth. Both effects are in 
accord with data for the solar neighbourhood. Finally, radial 
migration of stars will gradually flatten any radial metallic- 
ity gradient within the disc. We predict that there should 
be many stars with sub-solar metallicities at R ~ 5 kpc. 
The fact that metallicity gradients survive in disc galaxies 
in the teeth of such mixing puts a stronger requirement on 
the effectiveness of whatever process is responsible for their 
creation. 

Deflection of interstellar clouds by spiral arms at coro- 
tation should induce an anticyclonic circulation in the ISM 
that may already have been detected in observations of in- 
terstellar velocity fields. The circulation will bring into the 
solar neighbourhood gas from both low- and high-metallicity 
regions. Such importation of non-standard gas may account 
for the anomalously low metallicity of the Orion molecular 
cloud, and for the presence in the protosolar cloud of dust 
grains with large ratios n( 29 Si)/n( 28 Si) and n( 30 Si)/n( 28 Si), 
characteristic of highly nucleosynthetically processed mate- 
rial. 
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Our picture is incomplete, however, since we have im- 
plicitly assumed that the spiral behaviour manifested by the 
simulations correctly mimics the phenomenon in galaxies. 
Since the origin of spirals in galaxies, and in the simulations 
themselves, continues to lack a convincing explanation (e.g. 
Sellwood 2000) , this logical gap is likely to persist for some 
time. We regard the ability of simulations to predict the 
age- velocity dispersion relation of local disc stars (Carlberg 
& Sellwood 1985), the shape of the velocity ellipsoid (Jenk- 
ins & Binney 1990), the need for a dissipative component 
to maintain persistent spiral patterns (Sellwood & Carlberg 
1984), and the present success in accounting for the spread 
of metallicity of older stars, as indirect evidence to suggest 
that the recurrent transient patterns in the simulations do 
indeed mimic the behaviour in real galaxies. 
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APPENDIX: STAR-CLOUD SCATTERING 

Mass clumps on near-circular orbits induce a collective spi- 
ral response, or wake, from the surrounding disc (Julian & 
Toomre 1966). We regard such collective effects as part of 
the spiral structure of the disc which is the subject of the 
main part of the paper. Here we simply show that, when the 
collective wake is ignored, the interactions between stars and 
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the mass clump can be neglected as a possible source of ad- 
ditional angular momentum redistribution. 

For small radial oscillations, there is a convenient ex- 
pression for a star's epicycle energy in terms of the com- 
ponents ur and ii0 of its velocity with respect to the local 
standard of rest [BT eq. (7- 94a)]: 



Er = \{u R + 7 2 it 



4>)i 



(18) 



where 7 = 20./ k; i.e. 7 = v2 for a flat rotation curve. Small 
oscillations of a star perpendicular to the plane approxi- 
mately conserve the integral 



1„,2 
2 



(19) 



where $ is an effective gravitational potential. Following 
Spitzer & Schwarzschild (1953) we consider stars to be scat- 
tered in the impulse approximation, in which the scattering 
is complete in a distance small compared to both the epicycle 
amplitude and the smallest distance over which the Galac- 
tic potential changes significantly. In this approximation we 
have 



AE R = urAur + \{Au R f 



AE Z 



+-y 2 [u ct> Au lt> + ^(Au^) 2 ] 
: u z Au z + \(Au z )' 2 '■ 



(20) 



When a star encounters a massive object such as a GMC 
or a star cluster, the recoil of the object is negligible and it 
is natural to study the encounter in the frame of reference in 
which the scatterer is stationary throughout the encounter. 
If, moreover, the scatterer is moving on a circular orbit, 
this preferred reference frame is a steadily rotating one, and 
Jacobi's integral may be written [BT eq. (3-90)] 



E., 



IN 2 



$eff(r), 



(21) 



where v is the star's velocity in the cloud's rotating frame 
and $ e g is the sum of the gravitational and centrifugal po- 
tentials. Since we are treating encounters in the impulse ap- 
proximation, we may neglect the change in <& e ff between the 
in- and the out-states of the scattering event, and conclude 
from the invariance of Ej that the scattering merely redi- 
rects the velocity of the star in the rotating frame without 
changing its magnitude v. Moreover, at the radius of the 
cloud, u = v. We have therefore 

= Au = 2u r Aur + {Aur) 2 
+ 2u 4 ,Au 4> + (Au ) 2 + 2u z Au z + (Au z ) 2 . (22) 

Comparing this with equation ( pfj| ) and noting that AL = 
RAua, we obtain 



AL 



AE R + AE Z = ( 7 2 - 1) (u 4 + \AuA — 



R 



If we define 

Jrand = (Er 



E z 



then equation (|23|) can be written 
114, + |Au0 



AJr; 



(7 



-AL, 



(23) 



(24) 



(25) 



which is analogous to our result (^) for spiral waves. When 
a cloud scatters a star of the same specific angular momen- 
tum, u,f, — and the change in J ran d vanishes to first order 



because the interaction occurs at corotation for the star. 
Heating occurs when a cloud scatters a star that has a dif- 
ferent specific angular momentum because then uj, 7^ and 
the interaction occurs away from corotation. 

Since a small fraction of stars pass a cloud with — 0, 
most scattering by clouds is associated with heating. Hence 
we may use estimates of the amount of heating in discs 
to constrain the amount of angular momentum transfer for 
which clouds are responsible. Equation »25\) implies that 



AJ r 

and ^ 

Vc ' 



AL 



(26) 



Comparing this result with equation (|6j) for 
W0<5O/7km s _1 , we see that star-cloud encounters 
are about three times more effective at changing L for a 
given increment in random velocity than two-arm spirals 
at Lindblad resonance. Jenkins & Binney (1990) con- 
cluded that ((AJi ? ) 2 ) wavc ~ 90((AJ fl ) 2 ) c i oud . Here angle 
brackets imply averages over all encounters for a star that 
has Er/E z > 1 and an epicycle radius that is not large 
compared with the typical radial wavenumber of spiral 
structure. We have that 

— ((A Jfj) 2 ) wavc ~ ((A Jj?) 2 )cloud ~ ((A Jrand) 2 } cloud 

J? 2 ( 7 2 - l) 2 



(27) 

x (m^Am^) 2 + m (Am^,) 3 + i(Au^) 4 ) c i oud . 

The ensemble average on the right hand side of this equation 
is dominated by the first term, since the second averages to 
something near zero and the third term cannot exceed j of 
the first and is likely to be much smaller. Since scattering 
events occur at all radial phases, we may use the approxi- 
mation (u 2 ^{Au^) 2 ) ~ (uj,) {(Aurj,) 2 } . With equation (g) and 
setting 7 = \/2, we have finally 



<(AL) 2 ) cloud 



90 m?{ul) 



{(AL) 2 } 

Lindblad ■ 



(28) 



Hence, ((AI) 2 ) dolld = ((AL) 2 ) L mdbiad when («4>} 1/2 = 
w c /(9.5m) ~ 12 km s _1 for v c — 220 km s" 1 and m — 2. 
Since 12 km s _1 is a fairly typical value of (u$) 1//2 for an old 
disc star, we conclude that the angular-momentum changes 
induced by clouds are comparable to those induced by waves 
at Lindblad resonance. Since we have shown that waves in- 
duce much larger changes in L near corotation, clouds are 
not significant drivers of radial migration overall. 
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